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ABSTRACT 

The radio-loud quasar SDSS J102623. 61+254259. 5, at a redshift z=5.3, is one of the most 
distant radio-loud objects. Since its radio flux exceeds 100 mJy at a few GHz, it is also one 
of the most powerful radio-loud sources. We propose that this source is a blazar, i.e. we are 
seeing its jet at a small viewing angle. This claim is based on the spectral energy distribution 
of this source, and especially on its strong and hard X-ray spectrum, as seen by Swift, very 
typical of powerful blazars. Observations by the Gamma-Ray Burst Optical/Near-Infrared 
Detector (GROND) and by the Wide-field Infrared Survey Explorer (WISE) allow to establish 
the thermal nature of the emission in the near IR-optical band. Assuming that this is produced 
by a standard accretion disk, we derive that it emits a luminosity of ~ 9 x 10 46 erg s _1 
and that the black hole has a mass between 2 and 5 billion solar masses. This poses interesting 
constraints on the mass function of heavy (> 10 9 M Q ) black holes at high redshifts. 
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1 INTRODUCTION 

Blazars are Active Galactic Nuclei (AGN) whose relativistic jets 
are seen at a small angle from our line of sight, smaller than 1 /T, 
where T is the bulk Lorentz factor of the emitting plasma. The ra- 
diation emitted by their jets is therefore strongly boosted, making 
them visible at high redshifts. The most distant blazar known is 
Q0906+6930 (Romani et al. 2004; Romani 2006) with a redshift 
z = 5.47. The spectral energy distribution (SED) of this source 
reveals both the thermal (i.e. strong optical emission lines and con- 
tinuum) and boosted non-thermal (dominating in all bands but the 
near IR and optical) components. This is a common characteristic 
of powerful blazars at high redshift: the non thermal emission is 
characterized by two broad humps, peaking at lower frequencies as 
the bolometric luminosity increases, making the underlying ther- 
mal emission visible in powerful blazars (since it stands between 
the two non-thermal humps; Ghisellini et al. 2010a; 2010b). The 
high energy hump, in these sources, peaks below 100 MeV, disfa- 
voring its detection by the Large Area Telescope (LAT) instrument 
onboard the Fermi satellite (Atwood et al. 2009). Instruments ob- 
serving in the hard X-ray band can instead have more chances to 
detect them, and in fact the Burst Alert Telescope (BAT) onboard 
the Swift satellite (Gehrels et al. 2004) have detected blazars up 
to larger redshifts than LAT (compare Ajello et al. 2009 and Ack- 
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ermann et al. 2011). Their X-ray spectrum is hard [i.e. ax 0.5, 
assuming F{v) oc v~ ax ], and this, together with a relatively strong 
X-ray to optical flux ratio, can be taken as a signature of the blazar 
nature of the source. 

In this letter we suggest that SDSS 1102623.61+254259.5, a 
radio-loud AGN at z = 5.3, is a blazar, i.e. the viewing angle is 
smaller than 1 /V. If so, then this source indicates the presence, at 
approximately the same redshift, of ~ 2r 2 other sources point- 
ing in other directions, having similar intrinsic physical properties, 
including the black hole mass. Note that the flatness of the radio 
spectrum is not sufficient to guarantee the blazar nature of SDSS 
1102623.61+254259.5, since slightly misaligned sources (with re- 
spect to 1/T) still show a flat radio spectrum, especially at rela- 
tively large radio frequencies (at z = 5.3, the observed 1.4 GHz 
flux corresponds to a rest frame frequency of ~9 GHz), where the 
radio lobe steep radio emission does not contribute. Evidences for 
its blazar nature are its very large radio-loudness and very large ra- 
dio luminosity. Observed by Swift, it revealed a strong and hard X- 
ray flux, confirming its blazar nature. Observations performed by 
the Gamma-Ray Burst Optical/Near-Infrared Detector (GROND; 
Greiner et al. 2008) and by the WISE satellite (Wright et al. 2010), 
together with the Sloan Digital Sky Survey (SDSS; York et al. 
2000) spectrum, allowed to constrain the properties of the thermal 
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Table 1. GROND AB observed magnitudes of B2 1023+25, taken 2012 April 16th (magnitudes not corrected for Galactic extinction). The first raw gives the 
effective wavelength of the filter. 



2 B2 1023+25 AS A BLAZAR CANDIDATE 

To obtain a sample of high redshift blazar candidates, we se- 
lected a list of quasars included in the SDSS DR7 Quasar cata- 
log (Schneider et al. 2010), that have been analysed by Shen et 
al. (2011), with z > 4 and radio-loudness R > 100, defined 
as R = F(5GHz)/.F(2500A) (rest frame) in the work by Shen 
et al. We required a high radio-loudness in order to preferentially 
select objects with the jet directed toward us. Through this selec- 
tion, we obtained a sample of 31 sources with z > 4 including 3 
with z > 5. The latter sources have no rest frame 2500 A flux and 
hence we calculated F(2500A) extrapolating from the continuum 
at A = 1350A. 

SDSS J102623.61+254259.5 (= B2 1023+25) is the best 
blazar candidate from this list, because of its very high redshift and 
extreme radio-loudness (R ~ 5200). The precise value of the red- 
shift is somewhat uncertain, since z — 5.284 in NED, taken from 
the 6th Data Release of the Sloan survey (DR6), while z — 5.3035 
in DR7, and z — 5.266 in DR8. In this paper we assumed therefore 
a value of z=5.3. This is a known quasar, observed for the first time 
in the B2 Catalog of radio sources (Colla et al. 1972). The radio po- 
sition of this source is RA=10h26m23.62s, DEC=+25d42m59.4s 
and, according to the CRATES All-Sky Survey (Healey et al. 
2007), the source has a flat spectrum with a 1 4 /4. 8 = 0.489 [be- 
tween 1.4 and 4.8 GHz; F(v) oc u~ a ] and a 1A/8A = 0.504 with 
a radio flux of Fi.4ghz = 260.7mJy and Fsaghz = 105.7 mjy. 
We retrieved a number of radio fluxes from archival datcQ. More- 
over, B2 1023+25 is included in the WISE All-Sky Source Cata- 
log with clear detections in the two bands at lower wavelengths 
of the instrument, i.e. A = 3.4/im and A = 4.6/im. 

We retrieved and analyzed data of FermifLKT (Atwood et al., 
2009) collected in the period 2008 August 4th - 2012 June 14th (~ 
3.8 yrs). Photons of class 2 ("source"), with energies in the range 
0.1-100 GeV, collected in a region of 10° radius from the source 
position, were analyzed by using LAT Science Tools v. 9.27.1, In- 
strument Response Function P7, and the corresponding background 
maps. The followed procedures are described in detail in the anal- 
ysis threads of the Fermi Science Support Centefj. Particularly, the 
sources B2 1023+25 and Q0906+693 were fitted with a power law 
model, together with any other source of the 2FGL catalog (Nolan 
et al., 2012) within 10° from their sky coordinates. The analysis 
resulted in no detection for both sources, with a 5a upper limit on 
the photon flux above 100 MeV of 3.0 x 10~ 9 ph cm" 2 s" 1 . 



1 ASI Science Data Center (ASDC): 

|http ://tools.asdc.asi.it/| 

~ Data retrieved from http : //irsa . ipac ■ caltech . edu/| 
3 http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 



2.1 GROND observations 

The 7-band GROND imager, mounted at the 2.2 m MPG/ESO 
telescope at La Silla Observatory (Chile), started observing 
B2 1023+25 on 2012, April 16th at 23:16:56 UTC. We car- 
ried out two 8-minutes observations simultaneously in all 7 
g' , r' , i', z' , J, H, K s bands for a total exposure time of 919 s in 
the optical as well as 960 s in the NIR bands. Observations were 
carried out at an average seeing of 1.2" and at an average airmass 
of 2. 1 . The source was clearly detected in all bands but g for which 
an upper limit of 24.0 (AB magnitude) can be given. 

The GROND optical and NIR image reduction and photome- 
try were performed using standard IRAF tasks (Tody 1993) simi- 
lar to the procedure described in Rruhler et al. (2008). A general 
model for the point-spread function (PSF) of each image was con- 
structed using bright field stars, and it was then fitted to the point 
source. The absolute calibration of the g' ,r' ,i' , z' bands was ob- 
tained with respect to the magnitudes of SDSS stars within the 
blazar field while the J, H, K s bands calibration was obtained with 
respect to magnitudes of the Two Micron All Sky Survey (2MASS) 
stars (Skrutskie et al. 2006). 

Tab.[T]reports the observed AB magnitudes, not corrected for 
the Galactic extinction of E(B - V) = 0.02 from Schlegel et al. 
(1998). 



2.2 Swift observations 

Since this is our best blazar candidate, we requested a ToO ob- 
servation from the Swift team, in order to have X-ray data to 
confirm our hypothesis. The Swift satellite observed B2 1023+25 
twice: on 2012, June 21st (ObsID: 00032500001) and on 2012, 
June 22nd (ObsID: 00032500002). Data of the X-ray Telescope 
(XRT, Burrows et al. 2005) and the Ultra Violet Optical Telescope 
(UVOT, Roming et al. 2005) were downloaded from HEASARC 
public archive, processed with the specific Swift software included 
in the package HEASo ft v. 6 . 1 2Q, and analysed. The calibra- 
tion database was updated on July 16, 2012. We did not consider 
the data of the Burst Alert Telescope (BAT, Barthelmy et al. 2005), 
given the weak X-ray flux. 

The total exposure on the XRT was of ~10 ks and resulted in 
26 counts. Given the low statistics, the fit with a power law model 
with Galactic absorption (7V H = 1.50 x 10 20 cm~ 2 , Kalberlaet al. 
2005) was done by using the likelihood (Cash 1979). The output 
parameters of the model were T — 1.1 ± 0.5 and an integrated 
observed flux F .3-iokcV = (1-7 ± 0.4) x 10~ 13 erg cm -2 s" 1 . 
The X-ray data displayed in the SED (Fig.O has been rebinned to 
have 3cr in each bin. 

UVOT observed the source with two filters only: U, with 8.4 



4 Including the XRT Data Analysis Software (XRTDAS) developed under 
the responsibility of the ASI Science Data Center (ASDC), Italy. 
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ks exposure, and UVW1, with 1.6 ks. Both resulted in no detection 
with 3<r upper limits of the observed magnitudes equal to 21.70 and 
20.76 mag, respectively. 



2.3 Estimate of the disc luminosity 

We can calculate the disc luminosity Ld in two ways. The first is 
by assuming that the IR-optical flux is produced by a standard, 
Shakura & Sunjaev (1973) accretion disc, and finding the best Ld 
that accounts for the data. In this respect, for B2 1023+25 we have 
a rather good set of IR-optical data showing that the peak of the 
disc emission is at frequencies just below the absorption caused 
by intervening clouds (see Fig.rj}. Since the overall Ld is about 
twice the vL v peak luminosity, we directly derive Ld =— 9 x 
10 46 ergs-\ 

The second method is through the broad emission lines, as- 
suming they re-emit a fraction C (i.e. their covering factor) of the 
ionizing luminosity, and assuming that the latter is provided by the 
accretion disc. The average value of C is around 0.1 (Baldwin & 
Netzer, 1978; Smith et al., 1981), with a rather large dispersion. B2 
1023+25 shows a prominent broad Lya, that unfortunately is ab- 
sorbed in its blue wing. Using only the red wing, and subtracting 
the continuum in the same frequency range (adopting a power law 
model), we derive L LyQ ~ 2 x 10 45 ergs -1 and a FWHM=4000 
km s~ . These are twice the values we measured for the red wing 
only, implicitly assuming that the line is intrinsically symmetric. 
The ratio between the entire BLR luminosity and L^ ya has been 
calculated by Francis et al. (1991) and Vanden Berk et al. (2001). 
The former gives a Z/BLR/LLya = 5.55, while the latter gives a 
L B LR./L Lya = 2.7. Therefore L BL r ~(5-10)x10 45 erg s" 1 . Fi- 
0.1, the estimate for the disc luminosity is 
1 . Although this estimate is affected by 



nally, assuming C 
L d ~(5-10)xl0 46 erg s 
relevant uncertainties (e.g. partially absorbed Lya line and uncer- 
tain value of C), the estimated range encompasses the value found 
with the first method. 



3 BLACK HOLE MASS ESTIMATE 

The flux and spectrum of the disc depend on the central black hole 
mass A/bh and on the accretion rate M, The latter is traced by the 
total disc luminosity Ld = r/Mc 2 , with the efficiency parameter 
r\ assumed to be r\ — 0.08. If we measure Ld (hence Af ), we are 
left with A/bh as the only remaining free parameter. For a given M, 
the black hole mass determines the peak frequency of the disc emis- 
sion. A larger mass implies a larger disc surface and hence a lower 
temperature emitting a given luminosity. Therefore, for a fixed Ld, 
a larger A/bh shifts the peak to lower frequencies. Then the best 
agreement with the data fixes A/bh- We adopt Ld = 9 x 10 46 erg 
s _1 , as discussed above, and calculate the overall spectrum from a 
standard Shakura & Sunjaev (1973) disc. 

Fig. shows the data in the IR-optical band (as labelled) and 
three disc emission spectra calculated assuming the same Ld and 
three different mass values: M BH = 4.5 x 1O 9 M ( green line), 
2.8 x 10 9 M Q (blue) and 1.8 x 1O 9 M (cyan). Fig.[I]shows that, 
outside this range of masses, the model does not satisfactory repro- 
duce the data any longer. 
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Figure 1. Optical-UV SED of B2 1023+25 in the rest frame, together with 
models of standard accretion disc emission. Data from WISE and GROND, 
and the SDSS spectrum are labelled. Green points are archival data taken 
from ASDC SED builder. The grey stripe indicates the uL u peak luminosity 
of the disc. We show the spectrum of three accretion disc models with the 
same luminosity and different Mbh: log(AfBH/A40)=9.65 (green), 9.45 
(blue) and 9.25 (cyan). Note that outside this range of masses, the model 
cannot fit satisfactory the data. 
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Figure 2. SED of B2 1023+25 together with the adopted model. Data 
from WISE, GROND, Swift/XRT and Fermi/LAT are labelled. Green points 
are archival data taken from ASDC SED builder. The solid green line 
is the synchrotron component of the model, the dashed black line is the 
torus+disc+X-ray corona component, the dot-dashed grey line is the Ex- 
ternal Compton contribution. The LAT upper limits is for 3.8 years, 5cr. 



© 2012 RAS, MNRAS 000,[T]f5] 



4 Sbarrato, T. 




Log v [Hz, obs] 



Figure 3. Comparison between the SED of B2 1023+25 (red symbols) and 
of Q0906+6930 (grey symbols). As can be seen, the two sources are very 
similar. The grey vertical stripe marks the [0.3-10 keV] energy range of 
XRT, while the yellow vertical stripe corresponds to the [5-80 keV] energy 
range of the nuStar satellite. The curved grey stripe corresponds to the sen- 
sitivity of FermiLAT after 1 yrs of operation (5o\ lower bound) and of 3 
months (10 cr, upper bounds). We show the sensitivity (3<r) of the NuStar 
satellite for an exposure of 1 Ms. 



4 OVERALL SPECTRAL ENERGY DISTRIBUTION 

Fig.[2]shows the overall SED of B2 1023+25, together with a fitting 
model. The latter is described in Ghisellini & Tavecchio (2009), 
and consists of one emitting zone, in which relativistic electrons, 
whose energy distribution is derived through a continuity equation, 
emit by the synchrotron and Inverse Compton processes. The ac- 
cretion disk component is accounted for, as well the infrared emis- 
sion reprocessed by a dusty torus and the X-ray emission produced 
by a hot thermal corona sandwiching the accretion disc. The pa- 
rameters adopted to fit the SED are very similar to other powerful 
blazars studied and interpreted with the same model (Ghisellini et 
al., 2010a; 2010b). The bulk Lorentz factor is F = 14, the size 
of the emitting region is R — 5 x 10 16 cm, located at a distance 
of 5 x 10 17 cm from a black hole of Mbh = 2.8 x 1O 9 M . The 
magnetic field is B — 2.1 G. The main emission process for the X- 
ray and 7-ray non-thermal continuum is the inverse Compton pro- 
cess off the emission line photons (i.e. External Compton). Since 
the emitting region is rather compact, its radio emission is self- 
absorbed (up to~600 GHz) and cannot account for the observed 
radio flux, that must necessarily come from much larger zones. 
The Poynting flux carried by the jet is Pb ~ 10 48 erg s _1 , and 
the power spent to produce the radiation we see is P r ~ 3 x 10 45 
erg s - . The power in the bulk motion of the emitting relativistic 
electrons is P c ~ 2 x 10 44 erg s _1 , while the power in cold protons 
(assuming one proton per relativistic electron) is P p ~ 5.5 x 10 46 
erg s _ . 



4.1 Comparison with Q0906+6930 

Fig.[3]shows how the SED of B2 1023+25 compares with the SED 
of Q0906+693. As can be seen, the two SEDs are very similar. 
The only, remarkable, difference could be in the 7-ray band, if 
the tentative EGRET association with Q0906+6930 is real. On the 
other hand, the upper limits of FennifLAT are the same for the two 
sources. This may suggest that the EGRET flux is not produced by 
Q0906+6930, although variations of even more than two orders of 
magnitude of the 7-ray flux are not uncommon for blazars (see e.g. 
Abdo et al. 201 1 ; Ghirlanda et al. 201 1). 

Fig.[3]shows, as vertical stripes, the energy bands of XRT and 
of NuStar (Harrison et al., 2010). We report also the pre-launch es- 
timated NuStar sensitivity for an exposure of 1 Ms. It is remarkable 
that both sources have an X-ray flux lying one order of magnitude 
above the NuStar sensitivity, making their detection possible even 
with a moderate exposure time, and even at the high energy limit of 
the instrument, i.e. at ~80 keV. Note that 80 keV correspond, in the 
rest frame of a z — 5.3 source, to ~500 keV, that could be rather 
close to the peak of the high energy emission. 



5 DISCUSSION AND CONCLUSIONS 

In this letter we propose that the radio-loud, high redshift quasar 
B2 1023+25 is a blazar. It is the second most distant blazar known 
up to now, with redshift z = 5.3, that corresponds to an age of the 
Universe of 10 9 yrs. 

Both its thermal and non-thermal components are very lumi- 
nous. In agreement with the blazar sequence (Fossati et al. 1998) 
the two broad non-thermal humps peak at small frequency. In par- 
ticular, the hard X-ray spectrum and the upper limit in the 7-ray 
band constrain the high energy component to peak in the MeV re- 
gion of the spectrum. Therefore this source, along with the similar 
other powerful blazars, should have a relatively large hard X-ray 
flux and would have been an ideal target for planned hard X-ray 
instruments, such as EXIST, Symbol-X and NHXM, and it is a good 
target to be observed with the only orbiting satellite with focussing 
hard X-ray telescopes, i.e. NuStar (Harrison et al., 2010). The great 
sensitivity over its energy range [5-80 keV] enables it to detect the 
hard X-ray spectrum of this source, even if it cannot directly ob- 
serve the peak of the high energy hump. 

Thanks to the good data coverage in the IR-optical range, 
given by the combination of the 7 simultaneous optical-near IR 
filters of GROND coupled with WISE and SDSS data, we were 
able to estimate the central black hole mass, that turned out to be 
between 2 and 5 billion solar masses. This relatively small range 
compares well with the virial estimates in general (that have at 
least a factor ~4 of uncertainties, see e.g. Vestergaard & Peter- 
son 2006), and, in particular, it is the only method applicable for 
this source, where only a partially absorbed Lya line is visible (and 
virial methods are not yet well calibrated in this case). Near-IR 
spectroscopy, and hence the detection of the broad CIV and Mgll 
lines, would help to improve the virial estimate of the black hole 
mass and be a valuable cross-check of our estimate, based on the 
disk-continuum. 

It is remarkable to find radio-loud AGN hosting BH that mas- 
sive at an age of the Universe of only 10 9 yrs. The discovery that 
B2 1023+25 is a blazar, together with that one of Q0906+693, in 
fact, implies the existence of an entire family of radio-loud, very 
massive quasars in the early Universe, whose jet is not aligned with 
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our line-of-sight. Moreover, it is possible to roughly estimate the 
number of these objects and their spatial density. 

The comoving density of heavy black holes at high redshifts 
of radio-loud sources has been studied by Volonteri et al. (2011), 
based on the 3 years BAT catalog and the blazar luminosity func- 
tion, in hard X-rays, derived by Ajello et al. (2009), and modified 
(beyond z = 4.3) by Ghisellini et al. (2010a). In the latter paper 
the observational constrain on the blazar density with black holes 
heavier than 1O 9 M0 in the redshift bin 5 < z < 6 was based on 
the detection of only one object: Q0906+6930. Assuming it was 
the only blazar in the entire sky in this redshift bin, Ghisellini et 
al. (2010a) derived a comoving density of 2.63 x 10~ 3 Gpc" 3 of 
blazars hosting an heavy black hole in the 5 < z < 6 bin (see Fig. 
15 in that paper). 

B2 1023+25 was selected in the SDSS catalog covered by 
FIRST observations, and the common area of the sky of these 
two surveys is 8770 square degrees. Its existence in this portion 
of the sky implies that the space density of blazars hosting a 
M BH > W 9 M G black hole is a factor 40000/8770 greater than 
estimated previously (increasing the corresponding comoving den- 
sity to 1.2 x 10~ 2 Gpc -3 ). In turn, the density of all radio-loud 
sources hosting a Mbh > W 9 Mq black hole must be a factor 2F 2 
larger, namely ~ 4.7(T/14) 2 Gpc -3 in the 5 < z < 6 bin. Note 
that this is a lower limit, since other high redshift blazars could 
exist, but not yet identified as such. 

Moreover, assuming that the ratio of radio-loud to radio-quiet 
AGN is ~ 0.1, the SMBHs with M BH > 1O 9 M have to be at 
least 10 4 just after 10 9 yrs from the birth of the Universe. 
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